Mon. Not. R. Astron. Soc. 000, [T]{6] (2010) Printed 30 September 2010 (MN M£X style file v2.2) 



The Impact of Stellar Model Spectra in Disk Detection 



J. A. Sinclair, Ch. Helling, J. S. Greaves 

SUPA, School of Physics and Astronomy, University of St. Andrews, North Haugh, St. Andrews, Fife, United Kingdom, KY16 9SS 



Accepted 2010 July ?. Received 2010 July ?; in original form 2010 July 



ABSTRACT 



o 

(N 

in 

(N 
P? 

Oh 

6 



> 

m 
oo 

ON 

o 
o 

> 

X 



We present a study of the impact of different model groups in the detection of cir- 
cumstellar debris disks. Almost all previous studies in this field have used KURUCZ 
(atlas9) model spectra to predict the stellar contribution to the flux at the wave- 
length of observation thus determining the existence of a disk excess. Only recently 
have other model groups or families like MARCS and NextGen (phoenix) become 
available to the same extent as ATLAS9. This study aims to determine whether the 
predicted stellar flux of a disk target can change with the choice of model family - can a 
disk excess be present in the use of one model family whilst being absent from another? 
A simple comparison of KURUCZ model spectra with MARCS and NextGen model 
spectra of identical stellar parameters was conducted and differences were present at 
near-infrared wavelengths. Model spectra often do not extend in wavelength to that of 
observation and therefore extrapolation of the spectrum is required. In extrapolation 
of model spectra to the Spitzer MIPS passbands, prediction of the stellar contribution 
differed by 5 % at 70 fim for F, G and early K spectral types with differences increas- 
ing to 15 % for early M dwarfs. Analysis of the Spitzer MIPS 24 //m observations of 
37 F, G and K solar-like stars in the Pleiades cluster was conducted. In using KURUCZ 
model spectra, 7 disk excesses were detected while only 3 and 4 excesses were detected 
in using marcs and NextGen (phoenix) model spectra respectively. 



1 INTRODUCTION 



The study of Protoplanetary and Debris Disks is a signifi- 
cant area of research today. It is uncertain whether the pro- 
cesses which formed our own solar system are typical of other 
planetary systems. Both protoplanetary and debris disks de- 
scribe planetary systems in the midst of formation/evolution 
and therefore are crucial to our understanding of planetary 
systems as a whole. In steady state, the mutual collisions 
of planetesimals replenish dust removed from the disk by 
radiation pressure and Poynting-Robertson drag. Planetesi- 
mals and regenerated dust around older stars describe debris 
disks in contrast to the first-generation dust around stars 
younger than 10 Myr old. 

Detection of a disk may be gained if infrared emission from 
an object is in excess of expectation from the stellar flux. 
This is determined numerically by fitting a model spectrum 
and extrapolating out to the wavelength of observation - 



largely unused and have only recently become readily avail- 
able. 

Standard ID model spectra are the result of a solution of 
the equation of radiative transfer, the hydrostatic equation, 
mixing-length theory and the equations describing the gas 
phase in chemical equilibrium. Ideally, every possible atomic 
and molecular transition should be included however this is 
computationally impractical and simplifications are neces- 
sary, marcs and NextGen model spectra are produced us- 



ing Opacity Sampling (Helling & J0rgensen ( 1998 1) where 



opacities are distributed statistically in wavelength giving 
rise to high resolution model spectra. KURUCZ model spectra 
however are produced using Opacity Distribution Functions 
(ODFs) where absorption cross-sections are binned produc- 



ing a smoother, lower resolution model spectrum ( Gustafs 



son et al. ( 1975 1, Ekberg et al. ( 1986 1). Discrepancies in line 



see the Appendix of Bryden et al. ( 2006 1 for details of this 



procedure. The excess ratio (F ba/Fp re d) is computed and 
detection of an infrared excess is concluded when this value 



exceeds the 3<r limit of the entire sample (Aumann et al 



(19841, Beichman et al. (2005 1). In almost all past and 



features are expected between model groups due to these 
different opacity treatments though such differences are not 
expected at the longer wavelengths at which disks are ob- 
served. The different model groups include opacities drawn 
from different literature. This is likely to introduce sizeable 
differences in the spectral appearance of the model groups; 
for a discussion, see J0rgensen (20031. Model spectra with 



present studies of debris disks, kurucz (ATLAS9) mod- 



els (Castelli & Kurucz (2004 1 ) are used for this procedure. 



Other model groups like marcs (jGustafsson et aLJ ([2008)) 
and NextGen (phoenix) (Hauschildt et al. (19991) remain 



identical stellar parameters but of different model families 
are compared in Section 2 of this paper to determine whether 
such differences exist. The model spectra often do not ex- 
tend to longer wavelengths at which observations of disks are 
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neighbourhood stars which typically exhibit metallicities of 
[Fe/H] = 0.25 (jValenti fc Fischerj ||2005|)). 



6000 

Effective Temperature (K) 



Figure 1. The available effective temperature and surface gravity 
combinations for KURUCZ (red triangles), MARCS (blue squares) 
and NextGen (green circles) spectra of solar metallicity. The 
ranges of effective temperature which approximately represent 
the A, F, G, K and M spectral types are also displayed. 



made therefore, fitting of a Rayleigh- Jeans tail and extrap- 
olation to the required wavelength is necessary (Section 3). 
NextGen models do extend beyond 100 /im however, the 
sparse wavelength distribution at these wavelengths make it 
necessary to fit a Rayleigh- Jeans tail for interpolation onto 
some filter bandpass. Section 4 applies the use of all three 
model families to 37 solar-like stars previously analysed for 



disk excesses in Sierchio et al. (2010) 



2 COMPARISON OF MODEL SPECTRA 

Spectra of all three model families were compiled. KURUC^] 
and NextGen^] model spectra with inclusion of a 2.0 km 
s _1 microturbulence, a convective mixing length, l/H — 1.25 
and no parameterization of convective overshoot were cho- 
sen (| Castelli fc Kurucz| p004| ), |Hauschildt et al.| fl999| )). 
marc^] model spectra with identical microturbulence but 
(the only available at the time) convective mixing, l/H — 
1.5 were chosen ( |Gustafsson et al.] (|2008)). KURUCZ and 
NextGen model groups provide accurate synthetic SEDs 
of solar-like stars though KURUCZ better reproduces early- 
type stars while NextGen is more accurate in the M regime 
(Bertone et al. (20041). The Marcs model group have also 



been shown to adequately match the stellar SEDs given their 
use in the photometry calibration of the Spitzer Space Tele- 
scope ( Decin et al.| |2004 )) . Fig. [I] displays the effective tem- 
perature and surface gravity combinations of all three model 
families. As shown, model spectra in all three families ade- 
quately span the F, G and K spectral types which dominate 
the stellar sample of disk surveys. Model spectra of solar 
metallicity are considered the most appropriate in the con- 
text of disk surveys using stellar samples of solar-like, solar 



1 http:/ /kurucz. harvard.edu/grids. html 

2 ftp:/ /phoenix. hs.uni-hamburg.de/NextGen/Spectra/ 
3 



A comparison of the model spectra across the model fami- 
lies was conducted. A kurucz model spectra of some effec- 
tive temperature, surface gravity (and solar metallicity) was 
compared with a marcs and NextGen model spectra with 
identical stellar parameters. The wavelength range and res- 
olution differed between model families therefore interpola- 
tion (Eq.[T]) was necessary. In comparing KURUCZ and MARCS 
spectra, the former were interpolated onto the wavelength 
distribution of the latter. For comparison with NextGen 
models, a few short wavelength points were omitted to al- 
low a common wavelength range and the KURUCZ spectra 
interpolated onto the NextGen wavelength distribution. A 
ratio of the flux |2| was computed. 

(Fi - F 4 _0 



Fi — Fi-i + (Aj — Ai 



(Aj — A;_i) 



Fj 

Fratio(A) = -=7 



(1) 



(2) 



http:/ /marcs. astro. uu.se/search.php 



where Ai, Fi represents the KURUCZ model spectrum while 
Aj, Fj represents the marcs or NextGen model spectrum 
where A_, is a wavelength intermediate to A^_i and A ; . 
fratio(A) was calculated for each set of synthetic spectra and 
subsequently averaged for all spectra of the F, G and K spec- 
tral types. 

Spectral fluxes of KURUCZ and MARCS spectra can differ by 
greater than an order of magnitude (Fig. |2j. In particular, 
the poorest agreement occurs at approximately 0.2 /mi. The 
trend of the flux ratio is dominated by delta-function-like 
features which likely are produced as a result of line features 
present in one model group but absent from the other. This 
is likely capturing the different opacity treatments. A higher 
level of agreement occurs between kurucz and NextGen 
model spectra (Fig. |2J. Line features remain present, again 
capturing the different opacity treatments. The trends of 
both comparison tend to unity with increasing wavelength 
showing better agreement of the synthetic spectra at longer 
wavelengths. The best agreement is apparent for the KURUCZ 
and NextGen models while there exist some fluctuation 
about unity for the kurucz/marcs comparison. The bump- 
like feature at approximately 10 /im in the Kurucz/Marcs 
comparison possibly arises from differences in H~ opacity, 
the dominant source at these wavelengths. The deviation 
from unity in the Kurucz/NextGen comparison arises from 
an unphysical anomaly present in the Kurucz T = 8000 
K, log(g) = 4.5, [M/H] = 0.0 model spectrum. This feature 
was masked in all subsequent analysis. These differences in 
the model spectra has obvious implications in stellar flux 
prediction and finding the excess ratio from which detection 
of a disk is decided. 



3 EXTRAPOLATION TO LONGER 
WAVELENGTHS 

The algebraic form of the Rayleigh- Jeans' flux (Frj(X,T) — 
2kTc/\ 4 ) did not provide a good fit of the long wavelength 
regions of the synthetic spectra. Frj(\, T e g) (where T e g is 
the effective temperature of the model spectrum) was offset 
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Figure 2. The ratio of the KURUCZ spectral flux with (a) MARCS flux and (b) NextGen flux against wavelength (pm) averaged for all 
log(g) = 4.5, solar metallicity model spectra of F, G and K spectral types. The scales are logarithmic and the dashed red line represents 
perfect agreement. 



from the synthetic spectrum by up to 40% at 5 /im. A sim- 
ple scaling of Frj to the synthetic spectrum did not offer 
a solution as the fall-off of the synthetic flux with wave- 
length deviated from A -4 . Such an offset and slope differ- 
ence remained when the full blackbody flux was instead fit. 
A straight line fit of the spectrum in logarithmic space was 
instead conducted. Eq. [3] was fit (allowing freedom of a and 
b which minimised x 2 ) and transformed back into normal 
space, giving a function of the form shown in Eq. [4] 



\og{F x ) = a log(A) + b 



10" 



(3) 



(4) 



where a describes the wavelength power index (approxi- 
mately 4) while b serves as a scale for the pattern. Significant 
differences in the parameters which minimised \ 2 were ap- 
parent - for example, a ~ 3.9, b ~ 6.66 provides a good fit 
for a KURUCZ T = 3500 K, log(g) = 4.5 and [M/H] = 0.0 
model spectrum while a « 4.03, b w 6.94 provides a good 
fit for a NextGen model spectrum with the same parame- 
ters. For each model spectrum, such a fit was performed, the 
spectrum extrapolated out to 24 /tin and 70 /im and the flux 
integrated over the corresponding Spitzer MIPS passbands 
which are sensitive to dust emission at 1 AU and outside 10 



AU from a solar-like star respectively (Ricke et al. ( 2004 1) . 
For each set of kurucz and marcs or NextGen model spec- 
tra identical in stellar parameters, the broadband fluxes of 
these passbands were compared. Fig. [3] displays a ratio of 
these fluxes for each set of synthetic spectra. Kurucz model 
spectra do not extend to the cooler effective temperatures 
of M dwarfs and therefore Marcs and NextGen model spec- 
tra are also compared in this way to determine the possible 
differences in flux prediction for such objects. 

At 24 /mi, the extrapolated fluxes of KURUCZ and MARCS 
model spectra exhibit good agreement (of approximately 2 
%) for F, G and early K spectral types. The trend deviates 
significantly from unity for cooler late K and early M spec- 
tral types with disagreement as high as 15 % in KURUCZ 



and marcs spectra of 3500 K. At 70 /im, the flux ratio 
of KURUCZ and MARCS model spectra is identical to that 
of 24 /tm though it has a larger deviation from unity (ap- 
proximately 5%). This suggests that the long wavelength re- 
gions of KURUCZ and MARCS synthetic spectra differ in slope 
and therefore their log-log straight line fits (Eq. [3| diverge 
with wavelength. In comparison, the kurucz/NextGen 
trend exhibits a higher level of agreement which is consis- 
tent with Fig. [2] At 24 /tm, the flux ratio almost perfectly 
centres on unity with large disagreement only apparent in 
the comparison of early M dwarfs. At 70 /im, a small de- 
viation from unity is apparent though agreement is still 
true within 3 % for F, G and K spectral types. Similarly, 
Marcs and NextGen model spectra of F, G and K spec- 
tral types exhibit good agreement (of approximately 4 %) at 
24 /im, deviating to approximately 8 % at 70 /im. Again, M 
dwarfs exhibit the worst agreement with differences as large 
as 20 %. The poor agreement of M dwarf model spectra 
likely arises from the challenges in determining the opacity- 
relevant effects of molecular species present in such atmo- 



spheres (Gustafsson ( 1997 1, J0rgensen (20031, see also sec- 



tion 2.2 of|Helling fc Lucas| ( 2009 1). Observational studies 



have noted discrepancies between the observed and synthetic 



colours of such objects (Gautier et al. (2007l). 



4 ANALYSIS OF 37 SOLAR-LIKE PLEIADES 
STARS 

The possible impact in the use of different model groups 
in flux prediction of disk observations has so far only been 
inferred from comparison of the model spectra. It is therefore 
ideal to physically apply using all three model spectra to a 
set of observations. 



A recent publication (Sierchio et al. (2010)) analysed the 



Spitzer MIPS 24 /im observations of approximately 70 stars 
in the Pleiades cluster. Using KURUCZ model spectra for pre- 
diction of the stellar flux at 24 /im, they concluded that 23 
of such stars harboured disks. Using a smaller subset of their 
target sample (due to limitations in the available photom- 
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Figure 3. A comparison of the Spitzer MIPS predicted fluxes against the effective temperature of the log(g) = 4.5 KURUCZ and 
marcs/NextGen model spectra compared. The KURUCz/MARCS comparison is shown as blue triangles, the Kurucz/NextGen comparison 
are represented by red circles and the Marcs/NextGen comparison shown as the green squares. There are comparably fewer points for 
the KURUCZ and NextGen comparison as fewer spectra had matching stellar parameters. The dashed horizontal lines mark agreement 
of the fluxes. 



etry) this study was reproduced using KURUCZ, MARCS and 
NextGen model spectra. 

Johnson B, V, R, I and 2MASS J, H, Ks photometry of each 
object was compiled where possible. Synthetic photometry 
of all spectra of all three model groups was searched and x 2 
calculated in each instance. The observed and synthetic pho- 
tometry differed significantly by the target distance mod- 
ulus. Some targets in the Pleiades sample had known ac- 
curate (Hipparcos) distances however, for consistency and 
to avoid further reducing the stellar sample of the stars 
to those only with parallax measurements, an alternative 
method was decided. A mean magnitude in the above bands 
was calculated for both the synthetic and observed photo- 
metric magnitudes. The difference in these mean values was 
considered to give a value of the distance modulus. U band 
photometry was highly offset from longer wavelength pho- 
tometry, as noted in the Appendix of |Bryden et al. (2006), 
and was therefore omitted in this procedure. Chromospheric 
activity was considered to be the most likely explanation for 
these offsets in the U band (T. Lloyd-Evans, A. C. Cameron, 
private communication, 2010). Having found the synthetic 
spectrum in all three model families which provided the best 
photometry fit, the 'distance modulus' was added to the 
synthetic Spitzer MIPS 24 /im magnitude to determine the 
predicted flux. 

For each object, the predicted fluxes using all three model 
families were compared. Fig. [4] displays a plot of the pre- 
dicted fluxes of the 37 objects using KURUCZ model spectra 
against that of marcs and NextGen models. As shown, 
the line of equality provides an adequate best-fit of each 
data set. Although many points lie very close to equality, 
there are also many outliers and therefore, flux predictions 
for some objects using different model families can be sig- 
nificantly different. 

For each object, the excess ratio (F b B / F pre d) was calculated 
forming three distributions of excess ratios for each model 



Comparison of Predicted Fluxes 



1 1 1 1 1 


/ 




A 


Kurucz ( 1 ) vs. 


V* 


Mores (2) / 


^ Kurucz ( 1 ) vs. 


NextGen (2) / 












/ 

m * • 




















/ 

/ 

/ 




i 


i 



2 4 6 

Predicted Flux 2 (mjy) 

Figure 4. A comparison of the predicted fluxes of the 37 tar- 
gets | |Sierchio et~aj~| ( |2010[ |) at 24 fim (in mjy) in using KURUCZ, 
MARCS and NextGen model spectra. Red circles show the com- 
parison for KURUCZ and MARCS models while blue triangles show 
that of KURUCZ and NextGen models. The dashed line represents 
equality and is not a fit of either data. set. 

family. These distributions were binned and the resulting 
histograms are displayed in Fig. [5] As shown, the excess 
ratios form a gaussian-like distribution with positive out- 
liers being candidate excesses. It is immediately apparent 
that these distributions are somewhat different despite rep- 
resenting the same sample of stars. Gaussian functions were 
fit to these distributions in order to gain a quantitative mea- 
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Figure 5. The distribution of excess ratios (F i, a / F prel {) at 
Spitzer MIPS 24 (im for the sample of 37 stars using KURUCZ 
(red), marcs (blue) and NextGen (green) model spectra as a 
means of predicting the stellar contribution. Gaussian fits are 
shown as solid black lines and the x-translation and deviation pa- 
rameters of the fit (p, and ct) are also displayed. The 3cr confidence 
intervals associated with the Gaussians are shown as the vertical 
dashed lines. AKII437 and HII1132 (with very significant excess 
ratios) are omitted for the purpose of sensible axes scales. 



sure of their width and translation. A narrow distribution 
with a translation of fi — 1.0 describes model spectra which 
accurately predict the flux of the target stars. While the 
MARCS distribution is best-centred on unity, it is the broad- 
est of the three distributions. KURUCZ model spectra form 
the narrowest distribution while the NextGen distribution 
has the largest offset from unity. These offsets have been 
previously noted of Pleiades stars - see Fig. 4 of IGorlova] 
et~aTl((2006l. 



Detection of an infrared excess is considered conclusive when 
the excess ratio of an object exceeds the upper 3<r confidence 



interval of the entire distribution ( Sierchio et al. ( 2010 1). Ex- 
cluding AKII 437 and HII 132 (which significantly outlie), 
it is found that 7 objects (including 2 borderline objects) in 
the KURUCZ distribution satisfy this condition while falling 



to only 3 and 4 objects of the marcs and NextGen distri- 
butions respectively. There is agreement in all model fami- 
lies that AKIA76 (T cff = 6310 K, determined from the B-V 
colour), HII 1766 (6720 K) and PELS 146 (5600 K) have 
disk excesses. AKII 383 (6200 K), PELS 20 (5600 K) and 
PELS 150 (6200 K) yield disk excesses in using KURUCZ 
model spectra while no excesses exist if marcs or NextGen 
spectra are used. In particular, the KURUCZ, MARCS and 
NextGen spectra whose synthetic photometry best fit the 
photometry of AKII 383 have identical stellar parameters 
(T = 6000 K, log(g) = 4.5 , M/H = [0.0]). 

5 SUMMARY 

Detection of a disk can be concluded when infrared flux 
is in excess of expectation from the star alone. Popularly, 
a model spectra is used to predict the stellar contribution 
and almost all previous studies in this field have used KU- 
RUCZ (ATLAS9) model spectra for this purpose. MARCS and 
NextGen model spectra however have remained largely un- 
used due to only recent availability. 

An initial comparison of solar metallicity KURUCZ, MARCS 
and NextGen model spectra of identical stellar parameters 
yielded differences as large as two orders of magnitude. In 
extrapolating these spectra out to longer wavelengths, pre- 
diction of Spitzer MIPS fluxes could differ by 5 % for F, 
G and K spectra types and as much as 15 % for early M 
dwarfs. 



A recent publication (Sierchio et al. ( 2010 1) studied the 
Spitzer MIPS 24 /im observations of 70 stars in the Pleiades 
cluster for possible disk excesses, using KURUCZ (ATLAS9) 
models for prediction of the stellar contribution. This analy- 
sis was initially reproduced for a subset of these stars and 7 
out of 37 were determined to have disk excesses. In repeat- 
ing this analysis using marcs and NextGen model spectra, 
only 3 and 4 disk excesses were detected respectively. The 
differences between the model groups are such that different 
disk detection rates are obtained. We encourage the use of 
more than one model family to determine the presence (or 
absence) of an infrared excess. 
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